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SUMMARY 


This paper describes a technique, which has been programmed 
\s a DMAP Alter to Rigid Format 3, for determining a stiffness 
natrix modification to obtain a specified eigenvalue for a 
structure. The stiffness matrix modifications allowable are 
:hose that can be described as the product of a single scalar 
variable and a matrix of constant coefficients input by the user. 
The program solves for the scalar variable multiplier which will 
field a specified eigenvalue for the complete structure (pro- 
vided it exists) , makes the modification to the stiffness matrix, 
and proceeds in Rigid Format 3 to obtain the eigenvalues and 
sigenvectors of the modified structure. 


INTRODUCTION 


The motivation for devising a technique for determining a 
Local stiffness modification to obtain a specified eigenvalue 
9teramed from several launch loads analyses performed at the 
Soddard Space Flight Center in which these analyses were updated 
jsing data from hardmount spacecraft vibration tests . Quite 
aften, spacecraft are attached to their launch vehicle via a 
Garmon type clamp band. Generally, the clamp band attaches 
the spacecraft to an adapter section which in turn is bolted to 
the launch vehicle. However, the stiffness of the clamp band is 
aften not known well enough to make an accurate analytical pre- 
diction of the fundamental mode of the spacecraft adapter struc- 
ture when cantilevered from the base of the adapter, as it is in 
the spacecraft vibration tests. Thus, the original launch loads 
analyses are updated to reflect these discrepancies once the 
modes of the spacecraft-adapter structure have been measured in 
tests . 
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Updating any finite element model to agree with modal data 
obtained from tests usually requires a trial and error process 
in which some local stiffness is adjusted until the fundamental 
mode of the model agrees with the test data. However, if a value 
for the local stiffness exists which will give the finite element 
model the measured eigenvalue, then this stiffness can be found 
analytically . 

The DMAP Alter presented computes the value of the stiffness 
(or stiffness change) and adds this to the original stiffness 
matrix for the finite element model. The program then proceeds 
in Rigid Format 3 to compute the remaining eigenvalues and 
eigenvectors for the finite element model. 


THEORETICAL DESCRIPTION 


In real eigenvalue analysis, NASTRAN solves for the eigen- 
values and eigenvectors for the analysis ,^or U g ) degrees of 
freedom from 


T K aa ' XM aaJ 'V “ 0 <D 

The stiffness matrix for the U a degrees of freedom is obtained 
from the original Ug degrees of freedom through the application 
of constraints and Guyan reduction. The stiffness matrix Kgg 
for the Ug degrees of freedom can be considered to be the sum 
of two matrices 




+ AK 


gg 


(2) 


where K__ contains the stiffnesses for the finite element model 

yy 0 

which will not be modified and AK contains all of those stiff- 
nesses that will be modified. The modification technique des- 
cribed in this paper is one in which the stiffnesses to be 
modified are all proportional to some scalar variable, which will 
be denoted as P . Thus, AK g g can be written as 


4 
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gg 


(3) 



PJK 


where K gg ar ® the values of the^Kgg coefficients per unit value 
of the scalar variable p. The K matrix could represent, for 
example, the portion of the f initeFelement model represented by 
several beam elements of the same cross section whose moment of 
inertia we wanted to vary. In this case, p would be the moment 
of inertia of those beams and k' would be the stiffness 
coefficients for these beams per y Qnit moment of inertia. 

In general, AKgg can be any portion of the finite element 
model whose stiffness coefficients vary proportionally to some 
known variable. This variable could not, therefore, be the 
thickness of plate elements since the bending stiffness varies 
as the cube of the thickness while the transverse shear and 
membrane stiffnesses vary with the first power of the thickness. 
If, however, the plates were pure bending plates (no membrane or 
transverse shear), then all of the stiffness coefficients would 
vary with the cube of the thickness and we would be able to 
express the stiffness of those plate elements by an equation of 
the type in equation (3) where 3 could be taken as the cube of 
the thickness or the bending rigidity D. 

Thus, considering only those applications in which the 
stiffness matrix for a portion of the structure can be represented 
as in equation (3) where p is a single scalar variable, the 
stiffness matrix for the complete structure (eq. (2)) becomes 


l 



(4) 


The stiffness matrix in equation (4) can be reduced to the 
analysis set of degrees of freedom U a through the application 
of multi and single point constraints and through the Guyan 
reduction of the omitted points as mentioned above. The only 
restriction in the DMAP Alter presented herein is that the 
degrees of freedom that have stiffnesses that will be modified 
are not allowed to belong to the "0“ set (omitted coordinates) . 


271 


I 


Following the normal procedures for reducing from the to 
the U degrees of freedom (with the restrictions mentioned afiove) , 
the eigenvalue problem as stated in equation (1) can be written 
as 



+ 0K 


aa 


XM ] 

aa J 


(U > 


= 0 


(5) 


The problem is to find a value of 0 that will result in one 
of the eigenvalues (usually the first nonzero eigenvalue) attain- 
ing a specified value, say X-^. Setting X equal to the specified 
value X 1 in equation (5) results in the equation 


r 

L 



l M aa 


) + 0K 



{U } = 0 

a 


( 6 ) 


In order for there to be a nontrivial solution to equation 
(6) , the determinant of the coefficient matrix must vanish , This 
will result in a polynomial in 0 equal to zero, that is. 


P(P) = 0 


Thus, the solution for the value of 3 that will provide a 
specified eigenvalue (provided such value of 0 exists) may be 
obtained by solving an eigenvalue problem, using equation (6) , 
for 0 , This can be readily accomplished in NASTRAN using the 
module READ by inputting to READ the matrix (K - X 1 M ) as 

d 9 ^ ■» 9 9 

the "stiffness" matrix and the matrix K* aa as the "mass" matrix. 
The resulting "eigenvalue" found by READ will be the value of 0 
that will provide the stiffness modification necessary for the 
structure to have the real eigenvalue X^. 

It should be pointed out that there is no guarantee that the 
process will always work. There may be no modification of the 
portion of the structure we are attempting to modify that will 
result in the specified eigenvalue X^. However, the analyst 
can often tell, by comparison of his original finite element 
modes with those obtained from tests, what portion of the model 
appears to be too stiff or too flexible. In these instances, 
the procedure outlined in this paper for determining the stiffness 
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Modification should relieve the analyst of the burden of making 
Arbitrary changes in the stiffnesses and solving repeated eigen- 
value problems until the model agrees with the test. Since the 
technique outlined is one in which a stiffness change is deter- 
mined which will provide one eigenvalue equal to a specified 
Value, it appears that it will be most useful when there is 
is agreement between the original model and test results in a 
fundamental mode. It should also be mentioned that the stiff- 
ness change, while providing a specified fundamental mode, will 
obviously yield higher modes different from those obtained from 
~he original or unmodified finite element model. There is no 
guarantee that these new higher modes will agree any better with 
! :he test modes than those from the original model. 

i 

S INPUT TO THE PROGRAM 


The data deck required to make a run to modify part of the 
Structure and obtain the resulting eigenvalues will be discussed 
|.n terms of changes to a normal deck for Rigid Format 3, real 
eigenvalue analysis. 


Case Control Deck 

Two subcases are required. In the first subcase, a METH0D 
:ard selects an EIGB bulk data card which will be used for the 
dgenvalue extraction for p . 

The second subcase contains the normal case control cards 
:hat the user would have in any Rigid Format 3 run including a 
LETH0D card which selects the EIGR bulk data card for the real 
tigenvalues X • The result of this subcase will be the normal 
eal eigenvt .ue analysis output with one of the modes equal to 
:he specified eigenvalue (to be specified in the Bulk Data Deck) . 


Bulk Data Deck 

. Input of the normal finite element model of the structure 
which would be used in a real eigenvalue analysis. From 
this finite element model the stiffness matrix K will be 

built by NASTRAN. This could be the identical car^s used to 
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describe the structure if an original modal analysis had 
been performed and the user were now rerunning it to modify 
part of the structure . In this case, the value of 0 deter- 
mined in the current run would be the change in stiffness 
of the modified part of the structure . Included in these 
cards, of coviree, is the EIGR card requested by subcase 2 
which will find all desired modes subsequent to the modi- 
fication. 

2. DMIG input of k' 

gg 

3. EIGB card requested by subcase 1 for finding the "eigenvalue" 
0 . The normalization for the eigenvector must be MASS . 

If the scalar variable multiplier of K gg is, for example, 
the moment of inertia of some of the beam elements, then 
the search range should be the range over which the user 
expects the change in this variable to lie (change with 
respect to the value that is in the finite element model 
in item 1) . 

4. A PARAM bulk data card with parameter name = FR3Q and value 
equal to the frequency (in Hz) of the mode the user wishes 
to specify. 


DMAP ALTER DESCRIPTION 


Appendix A lists the DMAP Alters to Rigid Format 3, Level 
15.1.0, required to solve for the stiffness modification, to as- 
semble the new stiffness matrix, and to proceed in Rigid Format 3 
to obtain all of the desired eigenvalues and eigenvectors of the 
modified system. Several of the Alter statements are discussed 
in the appendix to clarify their function. In general, all 
the DMAP modules used but on^ are standard DMAP modules de- 
scribed in the NASTRAN User's or Programmer's Manuals. The module 
SCALAR, however, is a new module written and added to NASTRAN at 
the Goddard Space Flight Center and will be an available DMAP 
module in level 16 when it is released. Basically, this is a 
module that accepts matrices as input and will output one 
coefficient of the matrix as a NASTRAN complex, single or double 
precision parameter that can be used, for example, in the DMAP 
module ADD to multiply other matrices by. This was needed since 
the only way the scalar value of 0 could be obtained as data that 
could be used in subsequent DMAP statements was in the matrix 


2Jb 



+ * 


%*. KHHK output from module GKAM following the eigenvalue extraction 
t . for 0 . The module SCALAR was used to extract p from KHHK. The 
er matrix KHHK is the "modal stiffness" matrix found from the 
El' eigenvalue run to obtain 0. If the normalization on the EIGB 
p* bulk data card requests normalization to unit modal mass, then 
Id- the coefficient in KHHK will be 0 . 


SAMPLE PROBLEMS 


L 


Using the DMAP Alter program, two sample problems have been 
run. Figure 1 shows a beam finite element model of the UK- 5 
spacecraft and adapter to be flown on the Scout vehicle. The 
spacecraft and adapter are attached via a Marmon clamp, which 
in this finite element model is modeled as a scalar spring. In 
the original analysis, the model contained no scalar spring 
element for the clamp band and the adapter and spacecraft were 
assumed rigidly connected. The fundamental bending mode 
obtained from this finite element model was 43 Hz. Subsequent 
tests of the system indicated that the first mode was at 33 Hz 
and that the Marmon clamp did not appear "infinitely" stiff. 
Thus, the model vas modified by including a spring between the 
adapter and spacecraft. The second run, made to determine the 
value that the spring should have to obtain a 33 Hz fjxst 
bending mode contained the following changes: 


removal of the MPC rigid constraint at the adapter/spacecraft 
interface that was used in the original analysis to Simula ~e 
zero bending flexibility at that joint 


addition of DMIG matrix input of a scalar spring stiffness 
matrix per unit value of stiffness: 

fO 0 ... 0 0 . . . 0 1 


0 0 . . . 1 -1 ... 0 

0 0 ... -1 1 ... 0 


0 0 


0 0 


where the nonzero values correspond to the sows and columns 
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represented by the grid points and rotational degrees of 
freedom to which the scalar spring connects 

3 . EIGB bulk data card to find k g (3 is k g in this problem) 
with eigenvector normalization to MASS 

4. FARAM FREQ bulk data card with value 33 Hz (complex single 
precision) 

The data deck for this run is listed in Appendix B. The 
output from subcase 1 gave the value of k g needed to obtain a 
33 Hz first bending mode, namely, 4.3 x 10^ N/m (24.5 X 10^ lb/in) . 
Subcase 2 then was executed to obtain the eigenvalues and eigen- 
vectors for the system with this spring in the model. The 
resulting eigenvalues were a 33 Hz first mode with the second 
mode changing, in this case, by only a few percent from that 
obtained from the original model. 

Figure 2 shows another problem run using the DMAP Alter. 

In this case, the structure is a stiffened plate simply supported 
on all four sides. The plate is stiffened with an I-beam whose 
area and offset distance are specified but whose moment of inertia 
(about the beam centroidal axis) may be varied. The problem is 
to determine the moment of inertia of the beam that will give a 
40 Hz first symmetric bending mode of the structure. The struc- 
ture was modeled with a 5x5 mesh of grid points equally spaced in 
one quadrant of the plate. The DMIG matrix K'gg in this problem 
consisted of the stiffness of the beams (due to the bending moment 
of inertia only) for all of the grid points to which the beams 
were attached. The Bulk Data input for the finite element model 
consisted of the normal input for such a structure but with zero 
bending inertia for the beams (the area and offset distance were 
input on the CBAR cards) . The first subcase solved for the 
moment of inertia of the beams that would result in a 40 Hz first 
symmetric bending mode of the structure. This was determined as 
855.8 cm^ (20.56 in 4 ) . Subcase 2 then proceeded to obtain the 
eigenvalues and eigenvectors of the modified system and it was 
determined that the first mode was at 40 Hz. 
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APPENDIX A 


DMAP ALTER FOR DETERMINING LOCAL STIFFNESS CHANGE 
TO OBTAIN A SPECIFIED EIGENVALUE 
(RIGID FORMAT 3) 


1 

ALTER 45 


2 

MTRXIN 

• MATPOOL*EQEX IN* SIL* /DKGGP* t/V *N» LUSET/V* 
C*N*Q/C*N*0 $ 

3 

SAVE 

NODKP $ 

4 

matgpr 

GPL*USET »SIL*DKGGP//C*N*G/ r *N*G S 

5 

ALTER 48 


6 

EQUIV 

DKGGP tDKNNP/MPCFl % 

7 

ALTER 58 


8 

MCE2 

USET*GM*DKGGP***/DKNNP*** S 

9 

ALTER 61 


10 

EQUIV 

DKNNP*DKFFP/S INGLE $ 

11 

ALTER 64 


12 

UPARTN 

USET*OKNNP/DKFFP* * */C*N*N/C*N*F/C*N*S $ 

13 

ALTER 67 


14 

EQUIV 

OKFFP* DKAAP/OMIT S 

15 

ALTER 70 


16 

UPARTN 

uset*dkffp/dkaap***/c*n*f/c*n*a/c»n*o s 

17 

ALTER 75i 

76 

18 

AOD 

MAA*/MAA1/C*Y*FREQ s 

19 

ADD 

MAA1*/MAA2/C*Y*FREQ $ 

20 

ADO 

MAA2*KAA/0AA/C*N. (39. 47842*0. 0)/C*N* (-1.0 

21 

OPO 

dynamics*gpl*sil*uset/gpld*silo*usetd. * * • 


iEEO* 

EGOYN/V*N*LUSET/V*N*LUSETD/V*N*NOfPL/V*N*MOOLT/ 
V*N*NOP$DL/V*N*NOFRL/V*N*NONLFT/V*N*NOTRL/ 
V*N*N0EED/C*N*123/V*N*N0UE $ 
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22 

SAVE 

NOEED $ 

23 

CONO 

ERk0R2tN0EED S 

24 

CHKPNT 

EED i 

25 

READ 

OAA*DKAAP* * *EEO*US£T*CASECC/LAMAK*PMl AK »MlK*Ofc IGSK/ 
C*N*BUCKLING/V*N*N£IGVK/C*N*1 % 

26 

SAVE 

NEIGV $ 

27 

OFP 

LAM/,K*OE IGSK* * * *//V»N*CARDNOK S 

28 

SAVE 

CARDNOK $ 

29 

GKAM 

• PHI AK » M I K • LAM AK » » • , * CASECC/MHHK . * KHHK » PH I DHK / 
C*N*-l/C*N*l/C»Y*LFREO=0,0/C»Y*HFREQ*0.0/C.N*-l/ 
C*N*-i/C*Nf-l/V»NtNOCUP/V»N*FMODE i 

30 

SCALAR 

KHHK//C*N*1/C*N*1/V»N*BF.TA $ 

31 

SAVE 

BETA S 

32 

ADD 

DKAAP*KAA/KAAT/V*N*BETA i 

33 

COND 

LBL6*R£ACT $ 

34 

RBMG1 

USET » K A AT * M A A/KLL * KLR * KRR * MLL * MLR ♦ MRR S 

35 

ALTER 

85*90 

36 

REAO 

kaat*maa*mr*dh*eed*uset *casecc/laha*phia*mi *oeigs/ 
C »N*M0DES/V*N*N£I6V/C*N* 2 S' 

37 

SAVE 

NEIGV S 

38 

CASE 

CASECC* /CASEX2/C*N*TR AN/V*Nt REPEAT T S 2/V»N* NOLOOP S 

39 

alter 

105*105 

40 

S0R2 

CASEX2 *CSTH* HPT * DI T* EQEXIN*S IL* * *BGPDP* LAMA* QG» 
PHIG*EST*?*OOGi*OPHlG*OESl*OEFl*PPHiG/C*N*REIG S 

41 

alter 

o 

<• 

o 

o 

42 

PLOT 

PLTPAR«GPSETS*ELSETS*CASEX2*BGPDT*E0EXIN*SIP» *PPHIG/ 
PLOTX2/V*N*NSIL/V»N*I.USET/V*N*JUHPPLOT/V*N*PlTFLG/ 
VtNiPFILE S 


43 
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DESCRIPTION OP DMAP ALTER STATEMENTS 


5-16. 


18-20 


25. 


MTRXIN reads DMIG cards which contain the coefficients 
of the K gg matrix input by the user. These are the 
stiffness coef ficients (for the portion of the structure 
which will be modified) per unit value of the parameter 
that they vary with. These can easily be determined by 
running Rigid Format 1, up through GP4, with the bulk 
data containing all grid points, coordinate systems, and 
elements for the portion of the model to be modified. 


These Alters perform the reduction on the K matrix 
at the same location in Rigid Format 3 that the reduc- 
tions are performed on the stiffness matrix for the 
remainder of the structure (K a _ ) . 

gg o 


Formulate the matrix K - X,M using the input parameter 

I aa l aa 

FREQ which is FREQ = — . That is, FREQ is the fre- 

quency in Hz of the mode we are specifying the eigenvalue 
for . 


Solve an eigenvalue problem for pi using the buckling 
option in READ. The resulting "modal stiffness" matrix, 
KHHK, which will be output from module GKAM, will contain 
P on the diagonal since the eigenvector normalization on 
the EIGB bulk data card is a normalization on unit modal 
mass . 


29. GKAM outputs the matrix KHHK. 

30. SCALAR (discussed above) extracts a value from KHHK and 
outputs it as a parameter (BETA) . 

32. ADD formulates the total stiffness K „ + PK* . 

aa^ aa 

36. READ extracts the eigenvalues and eigenvector of the 
modified system, one of which will be the specified 
eigenvalue X^. 


279 


■ i REPRODUCIBILITY OF THE ORIGINAL PAGE IS POOR 
W 


APPENDIX 8 


CASE CONTROL AND BULK DATA DECKS FOR UK-5 S/C - ADAPTER STIFFNESS MODIFICAT 


TITLE = UK5 SPACECRAFT AND EH SECTION 
SUBTITLE = CANTILE/ERED MODE SHAPES IL ATERAL) 

LABEL = STIFFNESS CALCULATION FOR CLAMP BAND FOR 33 HZ FIRST BENDING 
ECHO = UNSORT 
MPC = 52 
SUBCASE 1 
METHOD = 1 
SUBCASE 2 
METHOD = 2 
OUTPUT 

VECTOR = ALL 
ELFQRCE = ALL 
SPCF = ALL 
BEGIN BULK 

S 

$ LATERAL MODES 


S 

GRDSET 

BAROR 


0 . 


EIGB 

1 

INV 

S,+b 

5. *7 

1 

♦EIG1 

MASS 





EIGR 

2 

INV 

25. 

400. 

3 

♦EIG2 

max 





PARAM 

GRDPNT 

0 




PAR AM 

WTMASS 

•002591 




S EH SECTION 





GRID 

SOI 


47.77 

0. 

0. 

GRID 

502 


44. 

0. 

0. 

GRID 

503 


40. 

0. 

0. 

GRID 

504 


37.27 

0. 

0. 

CBAR 

5001 

5001 

502 

501 


CBAR 

5002 

5002 

503 

502 


CBAR 

5003 

5003 

504 

503 


PBAR 

5001 

5001 

2.634 

94.4 

94.4 

PBAR 

5002 

5001 

2.138 

59.5 

59.5 

PBAR 

5003 

5001 

1.710 

29.2 

29.2 

♦B5011 

♦B5021 

♦B5031 

♦B5012 

.185 

.185 




♦B5022 

.185 

• 185 




♦B5032 

• IBS 

.185 




MAT1 

5001 

l.*7 


• 3 



S CONSTRAIN S/C - ADAPTER INTERFACE GRID POINTS 
$ SAME EXCEPT IN ROTATIONAL DEGREE OF FREEDOM 

S 

MPC 


MPC 


51 

52 


504 

504 


1 

2 


1.0 

1*0 


601 

601 


1345 

I. 0. 

1 

3 


123456 


72.6 .989 
45.8 .989 
22.4 .989 


TO BE THE 


1 - 1.0 

2 - 1.0 


1 

1.-4 

1.-4 


♦EIG 

♦EIG 


♦ 850 

♦B50 

♦850 

♦850, 

♦850 

♦850 


t 
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$ UK5 SPACECRAFT 


*, . 

' v ‘ 




M **** 







{GRID 

601 


37.27 

0.0 0.0 



GRID 

6C2 


33.77 

0.0 



GRID 

603 


30.27 

0.0 



,grid 

604 


26.77 

0.0 



GRID 

605 


23.27 

0.0 



{GRID 

606 


19,77 

0.0 



GRID 

607 


16.27 

0.0 



iGRiD 

608 


12.77 

0.0 



iGRID 

609 


10.27 

0.0 



iCBAR 

6001 

6001 

602 

601 



C8AR 

6002 

6002 

603 

602 



CBAR 

6003 

6003 

604 

603 



CBAR 

6004 

6004 

60S 

604 



CBAR 

6005 

6005 

606 

605 



CBAR 

6006 

6006 

607 

606 



CBAR 

6007 

6007 

60S 

607 



iCBAR 

6008 

6008 

609 

608 



IPBAR 

6001 

6001 

3.062 

37.05 

7.297 

♦B60 1 1 

Ipbar 

6002 

6001 

3.441 

44.87 

4.072 

♦B6021 

IPBAR 

6003 

6001 

3.525 

64. 

8.861 

♦86031 

Ipbar 

6004 

6001 

3.22S 

89.1 

10.809 

♦B6041 

!PBAR 

6005 

6001 

3.075 

118.5 

13.492 

♦B6051 

Ipbar 

6006 

6001 

3.165 

156.5 

10.618 

♦86061 

iPBAR 

6007 

' 6001 

3.505 

199.5 

7.328 

♦86071 

Ipbar 

6008 

6001 

3.965 

238.6 

25.31 

♦B6081 


♦B6011 

♦B6021 

♦86031 

♦B604| 

♦B605I 

♦B6061 

♦B6071 

♦B6081 

♦86012 .1907 

♦B6022 .1697 

♦B6032 .1589 

♦B6042 .1591 

♦B60S2 .1515 

♦B6062 .1349 

.1207 
.1067 


♦86072 

♦B6082 

HAT1 

DMIG 

OHIO 

OHIO 

♦OKI 

♦0K2 

PARAM 

ENODATA 


6001 

OKGGP 

OKGGP 

OKGGP 

601 

601 

FREQ 


l.*7 

0 

504 

601 

6 

6 

33. 


6 

6 

6 

- 1.0 

1.0 

0. 


.3 

1 


2 

504 

S04 


♦B6012 

♦B6022 

♦B6032 

♦860*2 

♦B6052 

♦B6062 

♦B6072 

♦B6082 


6 

6 


1.0 

~ 1.0 


♦OKI 

♦0K2 
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GRID POINTS 



AOAPTER UK- 5 SPACECRAFT 



THEORETICAL MODE SHA 

FOR k t = 4.3XI0 9 N/M 

(24.5XIO® LB/IN 
f = 33HZ 

THEORETICAL MODF SHA 

FOR k, — ► 00 
f = 43HZ 

O MODE SHAPE FROM 
UK- 5 /ADAPTER VIBRATI 
TESTS 
f =33HZ 


Figure 1.- Clamp-band stiffness modification to obtain 33 Hz first 
bending mode for the UK- j spacecraft and adapter. 
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PLATE 



SIMPLY SUPPORTED 


ALL 4 SIOES 


A * A 


$.08 CM 
(2.0 IN.) 


i mmnmnm 


r- PLATE THICKNESS^ $.08 CM 

(2.0 IN.) 

T (STEEL) 


BEAM AREA = 8.87 CM 2 (STEEL) 
( 1.375 IN. 2 ) 


Figure 2 .- Stiff en*r El modification to obtain 40 Hz first 
mode for the simply supported stiffened plate. 
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